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ABSTRACT 

Aims. The present dynamical configuration of planets in binary star systems may not reflect their formation process since the binary 
orbit may have changed in the past after the planet formation process was completed. An observed binary system may have been part 
of a former hierarchical triple that became unstable after the planets completed their growth around the primary star. Alternatively, 
in a dense stellar environment even a single stellar encounter between the star pair and a singleton may singificantly alter the binary 
orbit. In both cases the planets we observe at present would have formed when the dynamical environment was different from the 
presently observed one. 

Methods. We have numerically integrated the trajectories of the stars (binary plus singleton) and of test planets to investigate the 
abovementioned mechanisms. The orbits of the planets have been computed to test the survival of a planetary system around the 
primary during the chaotic phase of the stars. 

Results. Our simulations show that the circumstellar environment during planetary formation around the primary was gravitationally 
less perturbed when the binary was part of a hierarchical triple because the binary was necessarely wider and, possibly, less eccentric. 
This circumstance has consequences for the planetary system in terms of orbital spacing, eccentricity, and mass of the individual 
planets. Even in the case of a single stellar encounter the present appearance of a planetary system in a binary may significantly differ 
from what it had while planet formation was ongoing. However, while in the case of instability of a triple the trend is always towards 
a tighter and more eccentric binary system, when a single stellar encounter affects the system the orbit of the binary can become wider 
and be circularized. 

Conclusions. Modelling the formation of a planetary system around a binary is a potentially complex task and an effort has to be 
made to look into its present dynamics for traces of a possible chaotic past. 

Key words. Planets and satellites: formation; Methods: N-body simulations; Celestial mechanics 

1. Introduction tificial viscosity the presence of a companion star may help to 

trigger the formation of giant planets by disk instability. 

j=j Among the many extrasolar planets discovered so far, more than xhe differe nt models for planetary formation are applied 

40 are gas giants detected in binary or even triple star systems. consider ing that the present orbital configuration of the sys- 

This is not an unexpected result since it is known that about tem sta r-planet(s)-star was unchanged. However, there are cases 

half of solar type stars near the sun have a companion and there where the orbit of me companion star may have been modified 

is evidence that the binary frequency might be higher among bet ween the time when planetary formation was ongoing and 

young stars ([Reipurth a nd Zinneckeij [19931 | Kohler and Leinertj me prese nt. If the binary separation was wider or the eccentric- 

U9981 | Ghez etaU U993J). Most of the binary systems hosting a ity loweri the p rot0 pi ane tary disk would have been a more quiet 

planet have a wide separation ranging from 100 to 1000 AU (see environment promoting the formation of planets. There are two 

| DesideraandBarbieri| | 2007| for a compilation of binary systems possible mechanisms that may have altered a past binary dynam- 

with planets). However, there are some systems like y Cephei, ical configuration into the presently observed one. 
Gl 86 and HD 41004A where the separation between the two 

stars is of the order of 20 AU or less. Planetary formation in „. , ... . 

. . . , , . . - c . . . , - We assume that the observed binary system was in the past 

these systems is perturbed to a significnat extent by the grav- „ .,, , • , , • , 

, ■ a £ . .,, • part of an unstable larger multiple system (triple or more), 

itational influence of the companion star possibly altering the f., .„ , ,. , • , • , , 

FFu^u — i t t i N onn/iK u u Zu * tx. c c We will concentrate here on hierarchical triples, the most 

process. Thebault et al. (2004) have shown that the formation of ... . ,. , , r ' 

f, , n x. • • JFTTi — i — r~fl I mni common multiple system, in which the inner binary is or- 

the planet in y Cephei via core-accretion (Pollack et al, 1996) , • , , ■ • , • , , ■ r™ , • ; , 

, m, . tx. ,. , — j bited by a third star on a wider orbit. The dynamical evolu- 

is unclear. The migration in the massive gaseous disk required . -\, . , . J . , . , 

. c , .... t t c »tt £ . , . . turn of these systems may be complex in particular in the 

to form a planet within 2-2.5 AU from the star leads to an orbit ■; , r , ,. 

, , , . . , , ^ , , presence of high inclinations between the binary and the 

much closer to the star than that observed. On the other hand, , ■ , 

F; — I J-irvrv^ k i. u ,u . • . j j i j i -.u singleton. We may have systems whose parameters change 

Boss ( 2006) has shown that in hydrodynamical models with ar- , •■■„,, ■ , 

nearly periodically by large amounts that are stable over long 

timescales. However, if the triple is dynamically unstable 

Send offprint requests to: F. Marzari and the system becomes chaotic, the most frequent outcome 
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(stellar collisions are rarer) is the ejection of a star leaving 
a binary system with smaller separation, because of orbital 
energy conservation, and higher eccentricity. The planet in a 
circumstellar orbit around the main star of the pair has a time 
window to form in a quiet disk before the onset of dynamical 
instability of the triple system and subsequent modification 
of the original binary orbit into a more eccentric and closer 
one. 

- If the binary system formed in a densely populated cluster it 
is likely that the star couple underwent an encounter with an- 
other star. During the close stellar passage the orbital param- 
eters of the binary system are permanently modified leading 
to a different configuration star-planet(s)-star. Even in this 
case the mean trend is towards more eccentric orbit s and, on 
average, to an orbita l shrinking of the binary pair (Heggie , 
119751 iLarsonl 120011) . 

Both these mechanims, besides altering the binary configu- 
ration, have also the potential to destabilize the circumprimary 
orbits of potential planets. However, as we will see in the fol- 
lowing, inner circumprimary orbits, like that of the planet in y 
Cephei, are little or not affected by the stellar encounters in most 
cases. 

In Sect. 2 of this paper we will describe how the decay of 
a hierarchical triple star system affects the orbital parameters of 
the binary system that was part of it and of a putative plane- 
tary system around the primary star of the system. To simplify 
the approach we will model only coplanar hierarchical systems 
even if inclined triples may be more common because of fre- 
quent encounters in star forming regions early on in stellar evo- 
lution. This restriction to coplanarity is suggested by the minor 
complexity of the star dynamics, compared to the high inclined 
systems, that allows us to ouline clearly the effects of the stellar 
instability on the planet environment. We will also investigate 
to what extent the close encounters between the secondary stars 
perturb a planetary system around the main star. In Sect. 3 we 
will analyse the effects of stellar encounters on a binary system 
with planets. We will present selected cases that illustrate the dif- 
ficulty of interpreting the formation of a planet in a binary sys- 
tem in light of its present orbital paramters. We will not perform 
a full exploration of the parameter space for the two mechanisms 
described in Sect. 2 and 3 which would be beyond of the scope 
of this paper. 



2. Jumping stars and planet survival 

According to dReipurthl l2000t iKrouoal 119951: |E arson 1 1200 11) 
most binaries and even single stars originate from the decay of 
multiple star systems. The wide range of orbital parameters ob- 
served among binaries might be a consequence of the chaotic 
dynamics of the primordial multiple systems. Dynamical inter- 
actions among the members of the system cause an exchange of 
energy and angular momentum often resulting in the ejection of 
a star from the system. This may occur in the early stages of 
the system evolution or it may take some time for instability to 
build up. We will consider a scenario in which a triple star sys- 
tem is stable long enough to allow planet formation around the 
primary star of the binary pair. The subsequent destabilization 
of the triple, because of the onset of chaotic behaviour, will end 
with the ejection of one star. The binary couple will have orbital 
parameters that are at present different from those of the pri- 
mordial binary in which planetary formation occurred. In some 
cases, during the stellar encounters an exchange can also occur 



between the original companion star and the outer one leading 
to a new mass ratio of the binary system. 

2.1. The numerical algorithm 

To analyse the dynamical evolution of triple star systems and 
of planetary orbits around the primary star we have assembled 
a numerical model of the hierarchical three-body system. We 
define two initial osculating Keplerian orbits for the triple, an 
inner one for the binary pair nt\ and mo and an outer one, de- 
fined in the center of mass of the pair, for the single star m^. 
Ten additional circular orbits for massless bodies starting at 1 
AU from the primary star and extending by equal steps up to 10 
AU are also computed. Their evolution allows us to evaluate to 
what extent the chaotic phase of stellar encounters affects a pu- 
tative planetary system around the primary star. The evolution of 
the system is c alculated with the numerical integrator RADAU 
(lEverhartl 1 19851) which properly handles close encounters be- 
tween the massive bodies. 

The parameter space of the hierarchical three-body problem 
is rather wide since in principle there are 12 degrees of freedom: 
the six parameters of the binary pair and the six ones of the sin- 
gleton. In this work we restrict ourself to a limited number of sig- 
nificant test cases showing how a detached binary system has its 
orbital elements changed when the triple is disintegrated and the 
single star escapes to infinity. A similar behaviour occurs also 
in non-coplanar systems when the dynamics becomes chaotic, 
however we will perform a full exploration of the inclined case 
in a forthcoming paper. We start with the stars in prograde and 
coplanar orbits with the semimajor axis of the binary a\, fixed to 
35 AU. The semimajor axis of the single star a s is instead regu- 
lar ly sampled starting with i n the e mpirical stability limit derived 
by lEggleton and Kiseleval d 1995b . We consider different values 
of eb, the eccentricity of the binary, and of e s , the eccentricity of 
the singleton. Because of the chaotic character of the dynamics, 
for each set of (at,, a s , eb, e s ) we perform 30 different simulations 
with random choices of the orbital angles. In this demonstrative 
set of simulations the masses of the stars are 1, 0.4 and 0.4 solar 
masses, respectively. 

2.2. The onset of instability 

In FigQ] we plot the timescale f,„ required for the onset of in- 
stability, and then close encounters between the single star and 
the binary, as a function of the semimajor axis of the singleton. 
In the first case (red triangles) the initial orbital eccentricities eb 
and e s are both set to 0.2. There is a trend towards higher val- 
ues of ti n for larger values of a s that can be roughly described 
by an exponential law. Since we stop the simulations after 50 
Myr because of too long integration times, the curve is biased at 
the end of the fit by vales that are all equal to 50 Myr. Despite 
this limitation of th e fit, at a first sight the analytical dynamical 
stabil ity threshold (Eggl eton and Kiseleval [i 995: Kis eleva et all 
11996!) . located at about a s — \50AU for the singleton, appears 
to be in good agreement with our numerical results if one looks 
at the fast growing trend of the instability time when a s is larger 
than 150 AU. In the second case shown in FigQ](blue triangles) 
we have incresed both the eccentricities eb and e s to 0.4. The in- 
stability onset occurs, as expected, faster than the previous case 
and larger values of a s are required to have long surviving sys- 
tems. To find systems whose f,„ is longer than 50 Myr we have 
to go bey ond a s ~ 250 AU. The stab i lity Umit set by the f or- 
mula of (Eggl eton and Kiseleval 1 19951 iKiseleva et all 1 19961) is 
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Fig. 1. Time of the first close encounter between the binary and 
the singleton vs. a s , the initial semimajor axis of the singleton. 
The red dots mark the cases with binary eccentricity e\, = 0.2 and 
singleton eccentricity e s = 0.2, the blue dots those with eb = 0.4 
and e s = 0.4. The continuous lines are least squares exponential 
fits to the data. 
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Fig. 2. Distribution of the binary orbital elements after the third 
star of the system has escaped on a hyperbolic orbit. The red dots 
label binaries originated from triples with eb = 0.2 and e s = 0.2, 
the green dots those from triples with eb — and e s = 0. The 
large dots in the plot represent the initial binary system. 



at about 235 AU and it seems to be slightly less accurate at high 
eccentricities of the stars. 

According to Boss (2006) disk instability can form gas 
giant planets in a few hundred years. In this scenario the 
time required by the gravitational perturbations among the 
stars to build up dynamical instability is not a crucial pa- 
rameter. When a s is beyond 70 AU the chaotic behaviour 
onset occurs after a few thousand years giving enough time 
for planets to form. Eventuall y, if the planetary system 
is ch a otic after its formation (' Weidenschilling and M arzaril 
1996t iMarzari and Weidenschillingl l2002t iMarzaril 12001 
Rasio and Ford, 19961) it might have enough time to evolve into a 
stable state before the star system becomes in turn chaotic. In the 
core accretion theory dPollack et all 119961) . the formation of gas 
giant planets requires a few million years. By inspecting Fig[T] 
we expect that only systems in the upper part of the plot are dy- 
namically stable long enough to grant a quiet environment for 
planets to grow around the primary star. However, even in some 
cases where the singleton is close to the binary, the dynamics 
turns chaotic only after some million years. 

In Figf2]we show the final orbital distribution in the (ab,eb) 
plane of the left-over binaries. While the initial values of eb and 
e s have a strong influence on the instability time ?,■„, they ap- 
pear non influential for the final orbital distribution of the bina- 
ries. The two different distributions shown in Figf2] correspond 
to ei, = 0.2 and e s = 0.2 (red dots) while the green dots mark 
systems where both the binary and the singleton are on circular 
orbits. There is no significant difference between the two dis- 
tributions that overlap in the (a-b, eb) plane. The final semimajor 
axis ab is substantially smaller than the initial value (abo - 35 
AU) as expected by orbital energy conservation as one of the 
stars escapes on a positive energy hyerbolic orbit. The spread- 
ing of a b depends both on the different initial values of a s and 
on the amount of energy taken away by the singleton. The fi- 
nal values of eb are randomly distributed and are concentrated at 
high eccentricities, in most cases larger than 0.4. The results of 
additional simulations with sampled values of e s and e;, show a 
similar distributions in the final orbital distribution of the binary 
star system. 



In Fig[3]we map the values of the semimajor axis a p of the 
outer planetary orbit that remains stable after the chaotic phase 
vs. the final value of cib in the case with eb = 0.2 and e s = 0.2. 
The color codes are linked to the number of cases that fall in 
each bin in a p (x-axis) and ab (y-axis). The repeated encoun- 
ters between the stars set a limit of about 4 AU from the primary 
star within which planetary orbits can survive. Beyond that value 
of a p , the gravitational disturbances produced during the chaotic 
evolution of the two outer stars destabilize any planetary body. 
It is remarkable that a significant number of binary systems with 
final values of ab around 19 AU and high eccentricity, preserve 
the planetary system, or part of it, that formed when the binary 
was part of a triple system with abo = 35 AU. y Cephei, Gl 
86 and HD 41004A might well belong to this class of systems. 
In 1% of our simulations even the innermost planetary orbit we 
considered (a p = 1 AU) is perturbed during the stellar chaotic 
phase and the planet is eventually ejected from the system. This 
implies that a fraction of planetary systems forming around bi- 
naries, part of a triple in the early stages of their evolution, may 
be totally destabilized during the stellar encounters. 

Our study is not exaustive in terms of exploration of the 
parameter space, however it gives a clear idea of the dynami- 
cal evolution of the system. It also suggests caution when mod- 
elling planetary formation in an observed binary system by sim- 
ply adopting its present orbital configuration. That configuration 
might indeed be the outcome of a complex and unpredictable 
chaotic evolution. 



3. Modification of star-planet-star configurations 
after a stellar encounter 

Within stellar clusters, close stellar encounters may disrupt bina- 
ries or, more fr equently, abruptly change their orbital elements 
dKroupal Q995). This mechanism is suspected to be one of the 
possible causes of the low frequency of binaries among low- 
mass field stars compared to that of youn g low-mass stars in 
star-forming associations dDucheneUl999l) . If a planet or a full 
planetary system formed around one of the stars of a binary, 
when a stellar encounter occurs the dynamical configuration of 
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Fig. 3. Histogram illustrating the number of systems in our sim- 
ulations that, at the end of the chaotic phase, fall into bins in a/,, 
the binary separation, and a p , the semimajor axis of the outer 
planetary orbit that survives the stellar chaotic phase. 

the system may be significantly changed dDalla Stellal [2005). 
We may see either of the following events: 

- The binary system is destroyed and from then on the sur- 
viving planets orbit a single star. In this case, the dynamical 
configuration we observe at present is not indicative of the 
formation process because of the large changes caused by 
the stellar encounter. 

- The binary system survives the encounter but its orbital pa- 
rameters and those of the planets are strongly altered. 

In both cases, any attempt to model the formation of an ob- 
served planetary system would face the problem of discriminat- 
ing which dynamical or physical features of the system are due 
to the formation process and which are related to the stellar en- 
counter. 

In this Section we investigate to what extent the dynamical 
environment of a planet orbiting a star in a binary system is af- 
fected by a stellar encounter that does not disrupt the binary. 
We adopt the same numerical model described in the previous 
section, but the singleton is now set on a hyperbolic orbit ap- 
proaching the binary on a plane tilted by 45° with respect to 
the binary orbital plane. We randomly sample the impact param- 
eter and eccentricity of the hyperbolic trajectory and we look 
at the orbital configuration of the final system star-planet-star 
when the encounter is over. In Fig|4]we show the orbital distri- 
bution in the (a,e) and (a,i) planes of the binary systems after the 
stellar encounter. A significant number of systems has a post- 
encounter semimajor axis significantly smaller than the initial 
one (abo = 35 AU) while during the decay of a triple system the 
final Aft cannot be lower than a value fixed by the conservation 
of the gravitational energy of the initially bounded system. 

The impact of a stellar encounter on planetary formation 
models is important: the planet(s) may have formed when the 
companion star was farther away and then less effective in per- 
turbing the p r otopla netary disk. For the core-accretion model 
(Pollacketaj|,[l996), this may be a crucial condition to allow 
planetesimal accumulation before the onset of strong binary per- 
turbations when the compan ion stars is injected on a closer orbit 
( The bault et al.L l2004l I2006I) . By inspecting the planetary orbits 
that are still bound to the primary star after the encounter shown 
in Fig(5j we observe a distribution very different from that of 
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Fig. 4. Distribution of the eccentricity and inclination of the bi- 
nary system after the hyperbolic stellar encounter. All the sys- 
tems shown in the plot allowed the survival of at least one planet 
around the primary star of the system. The large dot represents 
the initial binary configuration. 
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Fig. 5. Plot showing the semimajor axis a p of the outer planetary 
orbit that survives the stellar encounter vs. the post-encounter 
semimajor axis of the companion star ab. To be compared with 
FigS 



Figf5] We still have cases where planetary orbits up to 3 AU 
survive the encounter leaving the planets around a closer binary 
system than the initial one. When we look at y Cephei, Gl 86 and 
HD 41004A we may be observing a system that survived a close 
stellar encounter. However, it appears that the single encounter 
with a star on a hyperbolic orbit has less dramatic consequences 
on the stability of a circum-primary planetary system than the 
frequent close encounters between bounded stars that lead to the 
decay of a triple. Most of the planetary systems show orbits that 
are still stable up to 10 AU when the stellar encounter is over. 
Of course, at later times outer orbits far from the primary may 
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become unstabl e because of the long term p erturbations of the 
companion star dHolman and Wie gerB,ll999h . 

While the onset of instability in a triple system always leads 
to tighter and more eccentric binary systems, by inspecting Fig|4] 
we see systems where the orbit of the binary couple is broadened 
and, in a few cases, even circularized by the stellar encounter. 
In these cases we know that planetary formation was perturbed 
by the companion star in the initial stages, but the dynamical 
environment of the planetary system became significantly qui- 
eter after the stellar encounter. There are systems where the bi- 
nary companion is moved beyond some hundreds of AU from 
the main star. If the stellar encounter occurred after the planets 
formed, we might observe a broad binary system with a small 
planetary system around the main star. Trying to understand why 
the orbits of such a planetary system are eccentric or too close to 
the star without considering the possibility of a stellar encounter 
would lead to wrong conclusions. 

4. Conclusions 

Any attempt to model the formation and evolution of a planetary 
system around the component of a binary star system should ac- 
count for the possibility that the binary orbit has changed with 
time. A dynamical event like the destabilization of an original 
triple system or a close stellar encounter can significantly al- 
ter the binary orbit after the formation of planets around the 
primary star. A system like y Cephei where a planet orbits the 
main star at about 2 AU and is perturbed by a secondary star 
presently moving at about 20 AU might have had a different dy- 
namical configuration when the planet formed. Trying to inter- 
pret the present configuration of the system using the observed 
orbital and physical parameters may lead to misleading deduc- 
tions on the physical properties of the protoplanetary disk that 
generated the planets. Prior to the binary orbit change, the disk 
may have been more (or less) affected by the gravitational pull 
of the companion star that was on a different orbit. The indi- 
vidual masses and orbit separations of the planets that formed 
in this environment might not be well simulated by introducing 
in the planetary formation model the present parameters of the 
system. Planetesimal accretion, which is the first step of terres- 
trial planet formation and of the core-accretion model for giant 
planet growth may occur differently depending on the binary or- 
bit and companion star distance. Even gravitational instability 
might follow a different path if the companion star was farther 
away in the initial phase of the binary life. 

In a few cases a stellar encounter may even cause the strip- 
ping of the companion star from the binary. The outcome would 
be a single star with a planetary system that formed when the 
star was part of a binary system and then perturbed by the grav- 
ity field of the companion star. In addition, stellar encounters 
may even push the orbit of the companion star out of its origi- 
nal plane leading to a significant mutual inclination between the 
planet and the companion. This is an efficient mechanism to pro- 
duce planets in a Kozai resonance with the secondary star. 
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